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ABSTRACT 

Although the nature of the progenitor of Type la supernovae (SNe la) is still unclear, the 
single-degenerate (SD) channel for the progenitor is currently accepted, in which a carbon-oxygen 
white dwarf (CO WD) accretes hydrogen-rich material from its companion, increases its mass to 
the Chandrasekhar mass limit, and then explodes as a SN la. The companion may be a main 
sequence or a slightly evolved star (WD + MS), or a red-giant star (WD + RG). Incorporating 
the effect of mass-stripping and accretion-disk instability on the evolution of WD binary, we 
carried out binary stellar evolution calculations for more than 1600 close WD binaries. As a 
result, the initial parameter spaces for SNe la are presented in an orbital period-secondary mass 
(logPjjMj) plane. We confirmed that in a WD + MS system, the initial companion leading to 
SNe la may have mass from 1 M s to 5 M . The initial WD mass for SNe la from WD + MS 
channel is as low as 0.565 M®, while the lowest WD mass from the WD + RG channel is 1.0 M©. 
Adopting the results above, we studied the birth rate of SNe la via a binary population synthesis 
approach. We found that the Galactic SNe la birth rate from SD model is 2.55 — 2.9 x 10 -3 yr _1 
(including WD + He star channel), which is slightly smaller than that from observation. If a 
single starburst is assumed, the distribution of the delay time of SNe la from the SD model may 
be a weak bimodality, where WD + He channel contributes to SNe la with delay time shorter 
than 10 8 yr and WD + RG channel to those with age longer than 6 Gyr. 

Subject headings: binaries: close - stars: evolution - stars: dwarf novae - supernovae: general - white 
dwarfs 



1. INTRODUCTION 
1.1. Progenitor Models 

Type la supernovae (SNe la) play an impor- 
tant role in astrophysics, especially in cosmol- 
ogy. They appear to be good cosmological dis- 
tance indicators and are successfully applied in 
determining cosmological parameters (e.g. il and 
A), which leads to the discover y of the accelerat- 
ing expansion of the u niverse ( Riess et al. 19981: 
IPerlmutter et al.1 1 1999h . At present, SNe la are 
proposed to be cosmological probes for testing the 
evolution of the Dark Energy equation of state 
with time (jHowell et al.ll2009al) . However, the pro- 
genitor systems of SNe la have not yet been con- 
fidently identified ( Hillebrandt &: Niemever 2000t 
iLeibundgutl l2000h . although they show their im- 



portance in many astrophysical field, e.g. deter- 
mining cosmological parameters, studying galaxy 
evolution, understanding explosion mechanism of 
SNe la and constraining the the ory of binar y stel- 
lar evolution (see the review by lLiviolll999h . 

There is a consensus that SNe la result from 
the explosion of a carbon-oxygen whit e dwarf in 
a binary system (jHovle fc Fowler 1960h . The CO 
WD accretes material from its companion, in- 
creases mass to its maximum stable mass, and 
then explodes as a thermonuclear runaway. Al- 
most half of the WD mass is con verted into ra - 
dioactive nickel-56 in the explosion (Branch 2004), 
and the amount of nickel-56 determines th e max- 
imum luminosity of SNe la (jArnetti 1 1982T) . Ac- 
cording to the nature of the companions of the 
mass accreting white dwarfs, two basic scenar- 
ios for the progenitor of SN la have been dis- 
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cussed over the last three decades 
single degenerate (SD) 



One is the 



1973 



. _ , model dW helan fe Iben 
Nomoto. Thielemann fe Yokoil Il984l ). in 



which a CO WD increases its mass by accreting 
hydrogen- or helium-rich matter from its com- 
panion, and explodes when its mass approaches 
the Chandrasekhar mass limit. The companion 
may be a main-sequence or a slightly evolved star 
(WD+MS) or a red-giant s tar (WD+RG) or a he- 
lium star (WD + He star) ( Yungelson et al.l 1995 ; 
Li fe van den Heuvellll997l: lHachisu et al. Ill999allb[ 



Nomoto et al 



" Eaneer et al.l 2000: 



200, 

Han fc Podsiadlowskil " Effli boOfil: IChen. fc Lj 
20071 120091: iHanl 120081: iMeng. Chen fc Hanll2009l 
Meng fc Yangj 120091: iLii et al.l l2009t IWang et al 
2009©. An al ternative is the double degener- 



ate ( DP) model (|Iben fc Tutukovlll984l IWebbink 



1984), in which a system consisting of two CO 
WDs loses orbital angular momentum by gravi- 
tational wave radiation and merges finally. The 
merger may explode if the total mass of the sys- 
tem exceeds t he Chandrasekhar mass limit (see 
the reviews by Hillebrandt fc Niemever 2000 and 



lLeibundgutll2000l ). Although the channel is the 



oretically less favored, e.g. double WD mergers 
may lead to acc retion-induced collapses rathe r 
than to SNe la (jfflllebrandt fc Niemeveil |2000|) , 



it is premature to exclude the channel at present 



contribute to a few SNe la (Howell et al.l 


2006 


Branch! 2006 


; Quimbv. Hoflich & Wheeler 


2007 


Hicken et al.l 


2007; 


Yuan et al.l 


2007|). At present 



the single-degenerate Chandrasekhar model is 
widely accep ted, since it is supported by many 
observations ( Parthasarathv et al. 20071 ) 



1.2. Observations 

The SD model is supported by many observa- 
tions. For example, variable circumstellar absorp- 
tion lin es were observed in the spectra of SN la 
2006X (jPatat et al1l2007l) which indicates th e SD 



nature of its precursor. iPatat et al.l (|2007l ) sug- 



gested that the progenitor of SN 2006X is a WD + 
RG system based on the exp ansion velocity of the 
circumstella material, while lHachisu et al. (2008) 
argued a WD + MS nature for this SN la. Re- 
cently, two twins of SN 2006X were also reported 
(jBlondin et al.l 120091: ISimon et al1l2009l ). and the 
birth rate and age of the 2006X-like supernovae 
may be well explained by the WD + MS model 



( Mc ng. Yang fc Gena 12009a ) . IVoss fc Nelemans 



20081 ) suggested that SN 2007on is also possi- 
bly from a WD + MS channel. Moreover, sev- 
eral SD systems are possible progenitors of SN la. 
For example, supersoft X-ray sources (SSSs) were 
suggested as good candidates for the progenitors 
of SNe la ( Livio 1996:lKahabka fc van den Heuvell 



19971 : lHachisu fc Katoll2003allbl ). Some of the SSSs 
are WD + MS systems and some are WD + RG 



systems (|Di Stefano fc Kondl2003l) . A direct way 
to confirm the progenitor model is to search for 
the companion stars of SNe la in their remnants. 
The discovery of the potential companion of Ty- 
cho's supernova may h ave verified the reliability of 
the WD + MS model dRuiz-Lapuente et al.1l200 " 



Branch 200J; llhara et al] l2007t iHernandez et al.l 



20091) . 

Measuring the delay time of SNe la (DD, be- 
tween the episode of star formation producing pro- 
genitor systems and the occurrence of SNe la) is 
a very important way to constrain the progeni- 
tor system. Although t he DD of most of SNe 
la is between 0.3-2 Gyr (|Schawinskill2009l ). there 
are still SNe la with long delay times (older then 
10 Gyr inferred from SNe la in elliptical galax - 
ies in the local universe, iMannucci et al.l 12005 ) 



or extre mely short delay times (shorter than 0.1- 
0.3 Gy r IMannucci et al.l 120061: Ischawinskil 120091: 
Raskin! 120091 ). Some observational results, i.e. 
the strong enhancement of the SN la birthrate 
in radio-loud early-type galaxies, the strong de- 
pendence of the SN la birthrate on the colors of 
the host galaxies, and t he evolution of the SN 
la birthrate with redshift ( Delia Valle et al. 2005t 
Mannucci et al. 20051 l2006h . may indicate a bi- 



modal distribution of delay time (DDT) , in which 
a part of the SNe la explode soon after starburst 
with a delay time less than 0.1-0.3 Gyr ('prompt' 
SNe la), while the rests have a much wider distri- 
bution with a delay time of ab out 3 Gyr ('tardy' 
SNe la IMannucci et~aT1 120061 ). 10% (weak bi- 
modality) to 50% (strong bimo dality) of all SN e 
la belong to the prompt SNe la ( Mannuccil 12008) . 
Whether this is due to the presence of two dif- 
ferent channels of explosion (such as single- and 
double-degenerate) or to a bimodal distribution of 
some parameter of the exploding systems (such as 
the mass ratio between the two stars of the binary 
system) is still unclear. 
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1.3. Situations 

Many works have concentrated on the SD 
model. To explain p rompt SNe la. variou s mod- 
els were designed. lHachisu et al.l (|2008l ) intro- 
duced a mass-stripping effect on a main-sequence 
(MS) or slightly evolved companion star by winds 



from a mass-accreting whit e dwarf (jHachisu et al 



Il996l) . IWang et al.l f|2009al ibh investigated WD 
He star channel which can explain SNe la with 
short delay times very well. iLii et al. ( 20091 ) even 
designed a symbiotic channel by assuming an as- 
pherical stellar wind with an equatorial disk. For 
the tardy SNe la, most of works focus on the 
WD + MS and WD + R G channel. Some au- 
thors (lHachisu et al.lll999allbl . 120081: iNomoto et ail 
Il999l [2OO3J) have "studied the SD channel by a 
simple analytical method to treat binary inter- 
actions. Such analytic prescriptions may not de- 
scribe some mass-transfer phases, es pecially those 



occur ri ng on a thermal time - scale (ILanger et al 



2000h . iLi fc van den Heuvell (jl997l ) studied the 



SD from detailed binary evolution calculation, 
while they considered two WD masses (1.0 and 
1.2 M Q ). ILanger et all (|2000l ) investigated the 
channel for metallicities Z — 0.001 and 0.02, but 
they only studied the case A evolution (mass 
transfer during core hydrog en burning phase). 
Han fc Podsiadlowski (20041) carried out a de- 



tailed study of this channel including case A and 
early case B (mass transfer occurs at Hertzsprung 
gap, HG) for Z = 02. Fo l lowing the study of 



Han fc Podsiadlowskil (|2004j ). iMeng. Chen fc Han 
(|2009l ) studied the WD + MS channel compre- 
hensively and systematically at various Z. Based 
on the works above, WD + MS channel is a very 
important channel to produce SNe la, but the 
channel may only account for about 1/3 SNe la 
observed and the de lay time for the channel is 
between 0.2-2 Gyr (lHan fc Podsiadlowskil l2004t 



Meng Chen fc Hani l2009h . Recently, |Xu_&Lj 



( 20091) proposed a WD + MS model with a ther- 
mal instable accretion disk, which can also pro- 
duce SNe la with long delay time. Although WD 
+ RG channel has been widely studied, it should 
be investigated carefully by detailed binary evo- 
lution calculations with the latest input physics 
since either the studies are based on a simple 
analytical method or focus on some special WD 
masses. In this paper, we want to construct a 
complete model, which includes WD + MS and 



WD + RG channels. 

In Section [21 we describe binary evolution 
model, and present the calculation results in Sec- 
tion |31 We describe our binary population syn- 
thesis (BPS) method in Section [4] and show the 
BPS results in Section [5j We give discussions in 
Section [6] and summarize our main conclusions in 
Section [71 

2. BINARY EVOLUTION CALCULA- 
TION 

2.1. Physical Input 



We use the stel lar evolution code of lEggleton 



(|197lL I1972L 119731 ) to calculate the binary evo- 
lutions of SD systems. The code has been up- 
dated with the l atest input physics over the last 
three decades (Han. Podsiadlowski fc Eggleton 



1993 IPols et al.llT995l 119981) . Roche lobe over 
flow (RLO F) is trea t ed w ithin the code de- 
scribed by lHan et al.l f)2000h . We set the ratio 
of mixing length to local pressure scale height, 
a = I /Hp, to 2.0, and set the conve ctive over- 
shoot i ng parameter Spy, to 0.12 dPols et al 



1997t ISchroder et al. 1997 ). which roughly cor- 
responds to an overshooting length of 0.25i?p. 
The solar metallicity is adopted here {Z = 0.02). 
Th e opacity table for th e met allicty is compiled 
bv IChen fc Toutl (120071) from llglesias fc Rogersl 



f|l996h and lAlexander fc Ferguson! (|l994l ). 



2.2. Accretion Disk &; Accretion Rates 

In an initial binary system, the companion fills 
its Roche lobe at MS or during HG or RG, and 
mass transfer occurs. If the mass transfer is dy- 
namically stable, the transferred material forms 
a disk surrounding the WD. If the effective tem- 
perature in the disk is higher than the hydrogen 
ionization temperature (~ 6500 K), the disk is 
thermally stable, otherwise an unstable disk is 
expected (|Osakilll996l : Ivan Paraduslll996l ). This 
corresponds to a critical mass-transfer rate, be- 
low which the disk becomes unstable. The critical 
mass-transfer rate depends on the orbital period 
of binary system: 

M c , th = 4.3 x IQ-^^^Mgyr- 1 , (1) 



wher e P nr h is the orbital p eriod unit in hour ()Osakil 
1996 : van Paradus 19961) . When mass-transfer 
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Fig. 1. — Three representative examples of binary evolution calculations. The solid, dashed, dash-dotted 
and dotted curves show the mass-transfer rate, M tr , the mass-growth rate of the CO WD, Mwd, the mass 
of the CO WD, Mwd, and stripping-off mass-lose rate from secondary, M str i p , respectively, in panels (b), 
(d) and (f). Dotted vertical lines in all panels and asterisks in panels (a), (c) and (e) represent the position 
where the WD is expected to explode as a SN la. 
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rate |M tr | exceeds the critical value, M c ,thi we as- 
sume the accretion disk is stable and the WD ac- 
cretes the transferred material smoothly at a rate 
M a = |M tr |. If |-Mtr| is lower than M Cjt h, mass is 
stored within the disk and no material is accreted 
by WD. When the materials stored in the disk 
reaches some critical amount, they are suddenly 
accreted onto the WD at a rate M a = \M tx \/d due 
to thermal instability , which may explain dwarf- 
nova outburst excellently (see the review bv lOsakil 
1996 for details about the disk-instability model). 
Here, d is duty cycle and is set to be 0.01 as did by 
Xu fc Lil (|2009h . Then, the accretion rate of WD 



(2) 



|M tr |, |M tr | > M Cjth , 
\M tI \/d, |M tr | < M c , th , 

where the typical timescale for the second accre- 
tion rate (in the case of the duty cycle) is equal to 
Mdisk/Ma, where Mdisk is the mass a ccumulated 
in the disk (see review by lLasotal < 200lh in details). 



2.3. WD Mass Growth 

Instead of solving stellar structu re equations of 



a WD, we ad opt the prescription o f lHachisu et al 
(|l999ah and lHachisu et all (|2008h on WDs ac 



creting hydrogen-rich material from their com- 
panions. The following is a brief introduction 
of this prescription. If WD accretion rate, M a , 
exceeds another critical value, M Cj h, we assume 
that a part of accreted hydrogen steadily burns 
on the surface of WD and is converted into he- 
lium at the rate of M Ci H- The unprocessed mat- 
ter is assumed to be lost from the system as 
an optically thick wind at a ra te of M W md = 
M a - M c ,h (|Hachisu et al.l Il996h . The optically 



thick wind may compact into the envelope of WD 
companion and strips off some hydrogen-rich ma- 
terial from the companion surface. The mass- 
stripping effect may attenuates the mass-transfer 
rate from the companion to the WD, and was suc- 
cessful ly adapted to explain some quasi-regular 
SSSs (lHachisu fc Katol l2003allbh . The mass-lose 
rate for the mass-stripping effect, M s trip is pro- 
portional to M w i nc i: 



AL 



strip 



ciM w 



(3) 



where c\ a constant. At present, the value of 
ci is very uncertain. To explain the position 
of a SSS (V Sge) in the initial orbital period- 
secondary mass (\0gP\M2) plane, c\ should be 



larger than (jHachisu et al.ll2008l ). lHachisu et al 
(|2008l ) checked the influence of c\ on the donor 
mass range and found that the range increases 
with c\. Considering primordial primary is more 
massive than primordial seco ndary, the c\ should 



be smaller 3 (see Fig. 3 in lHachisu et al.l 12008 
and the discussion in Section 15. In this pa- 
per, we assume rather arbitrarily that c\ = 1.5, 
which is the intermediate value between and 
3. Then, the mass-lose rate of the companion is 
M 2 = M t r ~ -MstrijB The material lost from the 
system may form circumstellar material (CSM), 
which is a possible origin of color excess of SNe la 
(|Meng et alJl2009h . 



The critical accretion rate is 



AL C h = 5.3x10 



-7( 1-7-*) 
X 



(AL 



WD 



-0.4)M Q yr-\ 
(4) 

where X is hydrogen mass fraction and Mwd is 
the mass of the accreting WD (mas s is in and 
mass-a ccretion rate is in M Q /yr, lHachisu et al 
1999ah . 



The following assumptions are adopted when 
M a is smaller than AL C h- (1) When M a is higher 
than iM c> H, the hydrogen- shell burning is steady 
and no mass is lost from the system. (2) When 
M a is lower than |M C ,H but higher than |M c ,h, 
a very weak shell flash is triggered but no mass is 
lost from the system. (3) When AL a is lower than 
|M q H, the hydrogen-shell flash is so strong that 
no material is accumulated on to the surface of the 
CO WE0. We define the growth rate of the mass of 
the helium layer under the hydrogen-burning shell 
as 

Mhc - f? H M a , (5) 

where rjn is the mass accumulation efficiency for 
hydrogen burning. According to the assumptions 
above, the values of 77H are: 




Ma > Mc,H, 

M C ,H > M a > |M CjH , 

M a < iM C)H . 



(6) 



1 Please note that M2 and Mtr are negative, while M s tiiv ls 
positive. 

2 Actually, when a WD undergos nova explosions, its mass 
is reduced. However, this should not affect our conclusions 
about SNe la. 
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Helium is ignited when a certain amount of he- 
lium is accumulated. If a He-flash occurs, some of 
the helium is blown off from the surface of the CO 
WD. Then, the mass growth rate of the CO WD, 
Mwd, is 



M WD = rj Hc M Hc = rfeemMa,, 



(7) 



where rjnc is the mass accumulation efficiency for 
helium-shell flashes, a nd its value is taken from 
Kato fc Hachisul <|2004h . 



We incorporated all the above prescription into 
Eggleton's stellar evolution code and followed the 
evolutions of both the mass donor and the ac- 
creting CO WD. The mass lost as optically thick 
wind is assumed to take away the specific orbital 
angular momentum of the accreting WD, while 
the stripped-off material by the optically thick 
wind is assumed to take away the specific orbital 
angular momentum of donor star. We calculated 
about 1600 binary systems, and obtained a large, 
dense model grid. M\ range from 0.6 M Q to 5 
M©; the initial masses of CO WDs, Af WD , from 

0. 565 Mq to 1.20 M ; the initial orbital periods of 
binary systems, P 1 , from the minimum value (at 
which a zero-age main-sequence (ZAMS) star fills 
its Roche lobe) to ~ 80 days. In the calculations, 
we assume that the WD explodes as a SN la when 
its mass reaches the Chandrasekhar mass limit, 

1. e. 1.378 M Q 

3. BINARY EVOLUTION RESULTS 

In Figs. [T]and[2] we present four representative 
examples and two extreme examples of our binary 
evolution calculations. The figures show the mass- 
transfer rate, M tr , the growth rate of CO WD, 
Mwd, the mass of CO WD, Mwd, stripping-off 
mass-lose rate, M str i p , the evolutionary track of 
donor star in the Hertzsprung-Russell (HR) dia- 
gram and the evolution of orbital period. 

3.1. Case Wind 

Panels (a) and (b) in Fig. 1 represent the evo- 
lution of a binary system with an initial mass of 
the donor star of M\ — 3.0M Q , an initial mass 
of the CO WD of M WD = 1.10M Q and an initial 
orbital period of log(P'/day) = 0.00. The donor 
star fills its Roche lobe at the MS stage which re- 
sults in a case A RLOF. The mass-transfer rate 



exceeds M c ,h soon after the onset of the RLOF. 
This results in a wind phase, where a part of the 
transferred mass is blown off as an optically thick 
wind, while the rest is accumulated on to the WD. 
The optically thick wind collides into the envelope 
of the companion and strips off some hydrogen- 
rich material from the surface of the companion. 
When the mass reaches M^, = 1.378Af Q , where 
the WD is assumed to explode as a SN la, the 
system is still in the wind phase (see the thick 
dotted line). At this point, the mass of the com- 
panion is M 2 SN = 1.9521M and the orbital pe- 
riod is log(P'/day) = —0.2050. Since the mass- 
transfer rate from the secondary continuously ex- 
ceeds M Ci h until the WD explode s as an SN la, 



we ca ll this "Case Wind" as did i n lHachisu et aL. 
(l2008h and iMeng. Yang fc Gend (l2009bh . Case 



Wind is realized in the region of 2.3M G 
3.8M Q and P l < 5 d ays for Ml 



WD 



< M\ < 
1.0 and 



l.ltffl. As noticed b v lHachisu et all (|2008l ) and 
Meng. Yang fc Gend (|2009bh . no Case Wind ex- 



ists for M WD < 0.9M Q (see also Figs. [3] and [p. 

Before supernova explosion, Case Wind may 
be observe d as quasi-r egular transient SSS such 
as V Sge ( Katol l2009'). The hydrogen- rich ma- 



terial lost from the system may form CSM and 
contribute to the intrinsic color excess of SNe la 
(jMeng et al. 2009). The material loses from the 



system via two ways: as optically thick wind from 
the surface of WD or stripped-off material by the 
wind from the companion. The velocity of wind 
should be larger than the escape velocity of the 
WD ( an order of magnitude of 10 3 km s" 1 ), 
which is too large to contribute the intrinsic color 
excess of SNe la (|Meng et al.l l2009h . The veloc- 
ity of the stripp ed-off material is re latively low 
(~ 100 km s"\ lHachisu et all 120081 ). Since the 



WD explodes during the wind phase, the stripped- 
off material may form CSM very close the SN la 
and may be the possible orig in of supernoyae like 



2008; 



SN 2002ic and SN 2006X (|Hachisu et al 
Meng. Yang fc Gend [2009ah . It may be helpful 
for explaining the pr operties of SNe la remnants 



(|Badenes et al.ll2007h . 



3.2. Case Calm 

Panels (c) and (d) in Fig. 1 show another ex- 
ample for an initial system with M\ = 2.50M Q , 
M WD = 0.90M Q , log(PVday) = 0.20. Similar 
to the previous example, the companion also fills 
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its Roche lobe at the MS stage and the system 
experiences a wind phase after the onset of the 
RLOF. After mass ratio reverse, |M tr drops un- 
til below M Cj H) but is still higher than |M Ci h, i.e. 
the optically thick wind stops and hydrogen-shell 
burning is stable. During this phase, WD reaches 
= 1.378M©, where M| N = 1.3773M© and 
log(P'/day) = —0.0460. The main difference be- 
tween this example and Case Wind is that at 
the time of the explosion, the system is in the 
stable hydrogen-burning phase after the optically 
thick wind phase while Case Wind is still in wind 
phase. Since WD undergoes steady hydrogen 
burning at the time of SN la explosion, we cal l 
this "Case Calm" as did in lHachisu et al.l (j2008l ) 
and lMeng. Yang fc Gengl (|2009bl ). Case Calm can 



be realized when > O.8OM (see also Figs. 

E] and HI). 

Before supernova explosion, C ase Calm may 
be observed as persistent SSS (jHachisu et al. 
20081: iKatol [2009). The materials lost as wind 



or stripped off by the wind form CSM, but they 
have been dispersed too far to be detected imme- 
diately after SN la explosion. We even can not 
expect radio or X-ray emission until at lease 10 
- 100 yr after the explosion. Then, Case calm 
would show the properties of "n o rmal" SNe la 



(Bra nch. Fisher fc Nugent 
20081 ) . 



1993; Hachisu et al 



3.3. Case Nova 

The third example in panels (e) and (f) of Fig. 
1 represents an initial system with M\ = 2.50Mq, 
M^ D = 0.80M Q and log(Pyday) = 0.10. Its evo- 
lution is similar to that of Case Calm except that 
the WD explodes at the phase where \M C ^{ > 
|M t r| > |M c ,h- When M^ D = 1.378M Q , Mf N = 
l.O816M and log(Pyday) = -0.1809. Before su- 
pernova explosion, hydrogen-burning is unstable, 
and t he system may be observed as a recurrent 
nova (jHachisu et al.ll2008t lKatol l2009l) . Therefore 



we ca ll this "Case Nova" as did in | Hachisu et al.. 
(|2008h and iMeng. Yang fc Penal (l2009bl) . Case 



Nova can be realized when M^ D > OJOMgH 



3 Please note that the WD masses in Figs. [3] and [4] are the 
initial masses of WDs, while the cases of "Case Wind", 
"Case Calm" and "Case Nova" are defined according to the 
evolutional stage of a WD binary system at the moment of 
supernova explosion. 



The materials lost as wind or stripped off by 
the wind form CSM, but they have been dispersed 
too far to be detected immediately after the SN 
la explosion. It takes at lease 100 yr for super- 
nova ejecta to re ach the CSM. So, a "normal" SN 
la is expected feranch. Fisher fc Nugent 19931 : 



Hachisu et alJl2008h . 



3.4. Case DNova 

Panels (a) and (b) in Fig. 2 represent the evo- 
lution of a binary system with an initial mass of 
the donor star of M\ = 3.5M , an initial mass 
of the CO WD of M^pj = 1.10M S and an ini- 
tial orbital period of log(Pyday) = 0.50. Af- 
ter RLOF occurs, the system experiences winds 
phase soon. Then, the system lost much hydrogen- 
rich material. After mass-ratio reverse, the sys- 
tem experiences the wind, the stable and weakly 
unstable hydrogen-burning phase one by one, un- 
til |M t r| < |M c ,h and even |Af tr | < M c ,th- At 
this phase, M WD still dose not reach to 1.378M©. 
Then, the WD increases its mass gradually by ac- 
cretion from a thermal unstable disk. After about 
1.3 x 10 6 yr, M^U reaches to 1.378M™, where 



M| N = 



yr, lviy^Y) readies iu i.oioiviq, 
0.6150M Q and log(Pyday) = 0.7548. 



e, the s 


ystem may 


( Osaki 


1996). So, 



"Case DNova" . Case DNova is mainly realized 
when APjA^wd ~ 3 and the system is crossing 
HG at the onset of RLOF. For Case DNova, af- 
ter |M tr | < |Af c ,H, |M tr | may be still higher than 
Af Cjt h- During this phase, the hydrogen-burning 
is heavily unstable, and then, the sy stem shoul d 
show the properties of classical nova (<Katdl2009h . 



For the same reason to Case Calm and Case 
Nova, a "normal" SN la is expected for Case 
DNova. 

3.5. Case Instability 

Panels (e) and (f) in Fig. 2 illustrate a more 
extreme case where both the donor and the WD 
are relatively massive. The initial binary in this 



case are M^ JVD = 1.10Af Q , 



M\ = 4.3A/ Q and 



log(Pyday) = 0.30. After about 1 x 10 5 yr from 
the onset of RLOF, mass transfer becomes almost 
dynamically unstable and hence the mass-transfer 
rate increases sharply, only to drop once the mass 



ratio has been reversed. When M, 
the binary parameters are 



WD 



= 1.378M , 
0.9461M Q 



s 



and log(P 1 /day) = —0.0566. For a larger initial 
donor mass, e.g. M\ = 4.5M©, our calculations 
show that mass transfer is unstable, and such sys- 
tems experience a delayed dynamical instability 
( Hiellming fc Webbinllll987t (Podsiadlowski et all 
20021 ). So, we call this case "Case Instability" 
There are mainly two differences between Case In- 
stability and Case DNova. One is that Case Insta- 
bility occurs when mass transfer begin at MS stage 
while Case DNova during HG. The other is that 
for Case Instability, both the donor and the WD 
are relatively massive, while Case DNova has no 
such constraint. 

3.6. Case RGB 

Panels (c) and (d) in Fig. 2 show a case that 
RLOF begin at red giant branch (RGB). The pa- 
rameters of the initial system are = I.IOMq, 
M| = 0.8M Q andlog^/day) = 0.50. For the sys- 
tem, the mass transfer rate is always lower than 
M C) h- There is even a broken phase for mass trans- 
fer. The WD mainly accretes hydrogen-rich mate- 
rial from a thermal unstable disk. After about 
2.7 x 10 8 yr, M™ D reaches to 1.378M©, where 
M| N = O.3277M and log(PVday) = 1.2528. To 
distinguish with Case DNova and due to the RGB 
nature of the donor, we call this "Case RGB". 
Case RGB is only realized when M^ D > 1.00M o . 
In our results, all WD + RG systems explode at 
the phase of M tY < M c ,th, and then they should 
be observed as dwarf nova before SNe la. 

3.7. Initial Parameters for the Progenitor 
of SNe la 

To conveniently compare our results with pre- 
vious studies in literatures, we summarize the fi- 
nal outcomes of all the binary evolution calcula- 
tions in the initial orbital period-secondary mass 
(\0gP\M2) plane. Figs. [3] and 0] summarize 
the final outcome of our binary evolution calcu- 
lations in the (log P'jM^) plane. Filled symbols 
show the results leading to SNe la, where the 
shape of the symbols indicates the supernova ex- 
plosions in the optically thick wind phase (filled 
five-pointed star: |M tr | > M c ,h)> after the wind 
phase while in the stable hydrogen-burning phase 
(Filled squares: M c ,n > \M tI \ > |M c ,h) or m 
the unstable hydrogen-shell burning phase (Filled 
circles: §M C ,H > \M tI \ > JM c ,h)) and unstable 



disk phase (Filled triangles: |M tr < M Cit h and 
\M tI \/d > |M C: h). Systems experiencing nova 
explosions and never reaching the Chandrasekhar 
limit and systems experiencing a dynamical mass 
transfer are also indicated in the figures. 

In Figs. [3] - [5J we also present the contours 
of the initial parameters in which SNe la are ex- 
pected. In the figures, the left boundaries are de- 
termined by the radii of ZAMS stars, i.e. RLOF 
starts at zero age. the WD + MS systems be- 
yond the right boundaries will undergo dynam- 
ically unstable mass transfer at the base of the 
red giant branch (RGB), while WD + RG sys- 
tems beyond the right boundaries will experience 
strong hydrogen-shell flash. The upper boundaries 
are determined by the delayed dynamical instabil- 
ity or the strong hydrogen-shell flash. The lower 
boundaries are constrained by the condition that 
the WD mass accretion rate is larger than |M c ,h 
and that the secondaries have enough material to 
transfer on to CO WDs, which can then increase 
their masses to 1.378 M . 

In these figures, we can clearly see that WD + 
RG channel occurs only when > LOOM© 

and the region for WD + RG channel increases 
with initial WD mass and extends to longer pe- 
riod. However, the maximum initial period for 
WD + RG channel is ~ 25 days, because the mass- 
transfer rate is very high and most of transferred 
material loses from system as optically thick wind 
or strippe d-off material if the period is very large. 
Similar to Xu &: Lj| §009), the lower boundary of 
WD + MS channel extends to ~ 1M , because 
the effect of a thermally unstable disk on the evo- 
lution of WD binaries is included in this paper. 
The upper boundary for WD + MS channel move 
to ~ 5M , because we in clude the wi nd mass- 



stripping effect suggested bv lHachisu et al.1 (|2008f ) . 
The high donor mass results might contribute to 
young population of SNe la. 

In Fig. [51 we overlay the contours for SN la 
production in the (log P 1 , AQ) plane for initial WD 
masses of 0.60, 0.70, 0.75, 0.80, 0.90, 1.0, 1.1 and 
1.2 Mq. Note that the enclosed region almost van- 
ishes when M-yyj) = O.565M , which therefore sets 
the minimum WD mass for which SD model can 
produce a SN la. We wrote the contours leading to 
SNe la into a FORTRAN code, which can be used 
for population synthesis studies. One can contact 
X. Meng to request the code. 
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Fig. 3. — Final outcomes of binary evolution calculation in the initial orbital period-secondary mass 
(log P 1 , M\) plane, where P 1 is the initial orbital period and M\ is the initial mass of donor star (for different 
initial WD masses as indicated in each panel). Filled five-pointed stars indicate SN la explosions during an 
optically thick wind phase (|Mtr| > M c ,h)- Filled squares denote SN la explosions after the wind phase, 
where hydrogen-shell burning is stable (M Cj h > |Af tr | > Filled circles denote SN la explosions 

after the wind phase where hydrogen-shell burning is mildly unstable (^M c ,h > l-^trl > |-^c,h)- Filled 
triangles denote SN la explosions during unstable disk phase (|Mtr < -Mc.th and |M tr |/d > |M Cj h). Open 
circles indicate systems that experience novae explosion, preventing the CO WD from reaching 1.378 M 
(|Mtr| < A^c,th and \M tI \/d < |M Cj h), while crosses show the systems that are unstable to dynamical mass 
transfer. 
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4. BINARY POPULATION SYNTHESIS 

Adopting the results in Section[31 we have stud- 
ied the supernova frequency from the SD channel 
via the rapid binary evol ution code developed by 
Hurley et ail (|2000L |2002h . Here after, we use the 



word primordial to represent the binaries before 
the formation of SD systems and initial for SD 
systems. 

4.1. Common Envelope 

Common envelope (CE) is very important for 
the formation of a SD systems. We firstly in- 
troduce the treatment for CE in this paper. 
During binary evolution, the primordial mass 
ratio (primary to secondary) is crucial for the 
first mass transfer. If it is larger than a criti- 
cal mass ratio, q c , the first mass tran sfer is dy- 
nami cally unstable and a CE forms (jPaczynski 
Il976h . The ratio <? c varies with the evolution- 
ary state of the primordial primary at the onset 
of RLOF iHiellming fc Webbinkl Il987t IWebbink 



al.ll2002 



iPodsiadlowski et al 



2002; 



19881: [Han et ... _ _ . . . 
Chen fe Hanll2008l ). In this study, we adopt q, 
4.0 when the primary is on MS or crossing HG. 
This value is supported by detailed binary evolu- 
tion studies (|Han et alJkoOOt IChen fc Hani 12002 . 
2003). If the primordial primary is on FGB or 
AGB, we use 



[1.67-z + 2(^) 5 ]/2.13, 



(8) 



where is the core mass of primordial primary, 
and x = d In R f /d In Mf is the mass-radius expo- 
nent of primordial primary and varies with com- 
position. If the mass donors (primaries) are naked 
heli um giants, Qr = 0. 748 based on equation @ 



Hurley et al.ll2002l for details) 



(see 

Embedded in the CE are the dense core of the 
primordial primary and the primordial secondary. 
Due to frictional drag with the envelope, the orbit 
of the embedded binary decays and a large part of 
the orbital energy released i n the spiral-in proces s 
is injected into the envelope ( Livio fc Sokerl [~L988'). 
Here, we assume that the CE is ejected if 



a C EA£ rb > \E b 



ind 



(9) 



where Ai? or b is the orbital energy released, 
-E'bind is the binding energy of common enve- 
lope, and q;ce is CE ejection efficiency, i.e. 



the fraction of the released orbital energy used 
to eject the CE. Since the thermal energy in 
the envelope is not incorporated into the bind- 
ing energy, acE may be greater t han 1 (see 
Han. Podsiadlowski fc Eggleton 19951 for details 
about the thermal energy). In this paper, we set 
«ce to 1.0 or 3.0. 

4.2. Evolution Channels 

There are three channels to produce WD + MS 
systems according to the situation of the primary 
in a primordial system at the onset of the first 
RLOF. 

Case 1 (He star channel): the primordial pri- 
mary is in HG or on RGB at the onset of the 
first RLOF (i.e. case B evolution defined by 



Kippenhahn fc Weigerdll9671 ). In this case, a CE 
is formed because of a large mass ratio or a con- 
vective envelope of the mass donor. After the CE 
ejection (if it occurs), the mass donor becomes a 
helium star and continues to evolve. The helium 
star likely fills its Roche lobe again after the cen- 
tral helium is exhausted. Since the mass donor is 
much less massive than before, this RLOF is dy- 
namically s table, resulting in a c lose CO WD+MS 



namicany s table, r esulting m a close UU W-U+lvu 
system (see lNomoto et alJll999l 120031 for details) 



Case 2 (EAGB channel): the primordial pri- 
mary is in early asymptotic giant branch stage 
(EAGB) (i.e. helium is exhausted in the core, 
while thermal pulses have not yet started). A CE 
is formed because of dynamically unstable mass 
transfer. After the CE is ejected, the orbit decays 
and the primordial primary becomes a helium red 
giant (HeRG). The HeRG may fill its Roche lobe 
and start the second RLOF. Similar to the He star 
channel, this RLOF is stable and produces WD + 
MS systems after RLOF. 

Case 3 (TP AGB channel): the primordial pri- 
mary fills its Roche lobe at the thermal pulsing 
AGB (TPAGB) stage. Similar to the above two 
channels, a CE is formed during the RLOF. A CO 
WD + MS binary is produced after CE ejection. 

However, only one channel can form WD + 
RG systems and then produce SNe la, i.e. the 
TPAGB channel. After a WD + MS system is 
produced from the TPAGB, the companion con- 
tinues to evolve until RG stage. Then, a WD + 
RG system forms. 

The SD systems continue to evolve and the sec- 
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ondaries may also fill their Roche lobes at a stage 
and Roche lobe overflow (RLOF) starts. We as- 
sume that if the initial orbital period, P^, and 
the initial secondary mass, M^, of a SD system lo- 
cate in the appropriate regions in the (logP 1 , M^) 
plane for SNe la at the onset of RLOF, a SN la is 
then produced. 

4.3. Basic Parameters in Monte Carlo 
Simulation 

To investigate the birth rate of SNe la, we fol- 
lowed the evolution of 10 7 bi naries via Hurley's 
rapid binary evolution code (jHurlev et al. 2000, 
2002). The results of grid calculations in section 
[3] are incorporated into the code. The primordial 
binary samples are generated in a Monte Carlo 
way and a circular orbit is assumed for all bina- 
ries. The basic parameters for the simulations are 
as follows. 



(i) The initial mass function (IFM) of lMiller fe Scalo 
fll979h is adopted. The primordial primary is gen- 
erated according to the formula of lEggleton et al 
(Il989h 



M{ = 



0.19A 



(l-A)°- 75 + 0.032(1 -X] 



0.25 ' 



(10) 



where X is a random number in the range [0,1] and 
Mf is the mass of the primordial primary, which 
ranges from 0.1 M Q to 100 M Q . 

(ii) The mass ratio of the primordial com- 
ponents, q, is a very important parameter for 
binary evolution while its distribution is quite 
controversial. For simplicit y, we take a uni- 



form mass-ratio distribution ( Mazeh et al.l 11992; 



Goldberg fe Mazehlll994i) : 



n(q) = 1, 



<<?<!, (11) 



where q = M$/Mf. 

(iii) We assume that all stars are members of 
binary systems and that the distribution of sepa- 
rations is constant in log a for wide binaries and 
falls off smoothly at close separation: 



an(a) 



_ / a sep (a/a y 



a < ao; 
ao < a < oi, 



(12) 



where a sep 0.070, ao = 10i?Q, a\ — 5.75 x 
10 6 i? Q = 0.13pc and m w 1.2. This distribu- 
tion implies that the numbers of wide binary sys- 
tem per logarithmic interval are equal, and that 



approximately 50% of the stellar systems are bi- 
nary systems with orbital periods less than 100 yr 
(|Han. Podsiadlowski fc Eggleton|[l995l) . 

(iv)We simply assume a single starburst (i.e. 
1O 11 M0 in stars are produced one time) or a 
constant star formation rate S (SFR) over the 
last 15 Gyr. The constant SFR is calibrated so 
that one binary with Mi > 0.8M P) is born in 
the Galaxy each year (see Iben fc Tutukov 1984 



Han^ Podsiadlowski fc Eggletonll995 ; Hurley et al 
2002). From this calibration, we can get 5* = 5 



M yr" 1 (see also IWillems fc Kolbl 120041 ) . The 
constant star formation rate is consistent with the 



estimation of lTimmes et al.l (|1997l ) , which success- 
fully reproduces the 26 Al 1.809-MeV gamma-ray 
line and the core-collapse su pernova rate in the 
Galaxy ( Timmes et al.lll997 ). Actually, a galaxy 



would have a complicated star formation history, 
while the two choices here are extremes for sim- 
plicity. A constant SFR is a g ood approximation 



for spi ral galaxies like Galaxy ([Rocha-Pinto et al 



2000al lbT). while a single starburst is for elliptical 
galaxies. 

5. The RESULTS of BINARY POPULA- 
TION SYNTHESIS 

5.1. the Birth rates of SNe la 

Fig. [6] shows Galactic birth rates of SNe la for 
the SD channel. The simulations give a Galactic 
birth rate of 2.25 - 2.6 x 10~ 3 yr" 1 (thick lines), 

only slightly lower than the b irth rate inferred ob- 

serva t ionally (3-4 x 10~ 3 yr~ 1 , van den Bergh fc Tammann 



199ll ICappellaro fc Turattolll997h. If WD + He 



star channel is included (jWang et a l. 2009b), the 
Galactic birth rates of SNe la for the SD channel 
becomes 2.55 - 2.9 x 10~ 3 yr" 1 (thin lines). The 
WD + He star channel contributes to SNe la by 
about 10 per cent. 

Fig. [7] displays the evolution of birth rates of 
SNe la for a single starburst of 10 11 M . Most of 
the supernova explosions occur between 0.15 and 
2.5 Gyr after the burst, which may mean that S/S0 
galaxies with ages below 3 Gyr have a higher SNe 
la birth rate than th ose older galaxies (see obser- 
vatio nal results fromlGallagher et al.ll2008h . Simi- 
lar to lHan fc Podsiadlowskil (|2004l) , we also found 
that a high «ce leads to a systematically later ex- 
plosion time, because a high acE leads to wider 
WD binaries, and, as a consequence, it takes a 
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Fig. 5. — Regions in the initial orbital period- 
secondary mass plane (log P^M^) for WD bina- 
ries producing SNe la for initial WD masses of 
0.60, 0.70, 0.75, 0.80, 0.90, 1.00, 1.10 and 1.20 
Mq . Since the region almost vanishes for = 
0.565M Q , we do not show it in the figure. 




log(i/yr) 

Fig. 6. — The evolution of the birth rates of SNe 
la for a constant star formation rate (Z=0.02, 
SFR=5M Q yr _1 ). Solid and dashed lines show the 
cases with «ce = 1-0 and q;ce = 3.0, respectively. 
The thick lines are the results in this paper, while 
the thin lines are the results including th e WD + 
He star channel from Wang et al. ( 2009bl ). 



longer time for the secondary to evolve to fill its 
Roche lobe. Actually, because of the low binding 
energy of the common envelope and a long primor- 
dial orbital period, acE has a remarkable influence 
on CO + WD systems from the TPAGB channel. 
Generally, if a CE can be ejected, a low acE pro- 
duces a shorter orbital-period WD + MS system, 
which is more likely to fulfill the conditions for SNe 
la. Therefore, we see obvious contribution from 
the TPAGB channel when acE = 1-0, but no con- 
tribution from this chann el when «ce = 3.0 (see 



Meng. Chen &: Hanll2009l for details regarding the 



influence of «ce on the TPAGB channel). 

In Fig. [7j it is clear that there is a high delay- 
time tail extending to about 15 Gyr, since the WD 
+ RG system is also included in this paper. How- 
ever, the contribution of WD + RG channel to SNe 
la is much smaller than that of WD + MS chan- 
nel. The WD + RG channel mainly contributes 
to SNe la with delay time longer than 6 Gyr. 

(|2008l) claimed that 



Although IHachisu et al 



the model with mass-stripping effect may produce 
very young population of SNe la, there is no SNe 
la with delay time shorter than 10 8 yr in our BPS 
results. This result is determined by binary evolu- 
tion, i.e. primordial primary must be more mas- 
sive than primordial secondary. Let us consider 



a simple example. For a CO WD o f M, 



M(7) , i ts progenitor mass is ~ 6 M (jUmeda et al 



wn 



1.1 



1999; iMeng et al.ll2008l ). Then, the upper limit 



of primordial secondary mass for an initial CO 



WD binary system with M^n — 1.1 Mq is ~ 6 
Mq . When primordial primary evolves to TPAGB 
phase, a CE forms because of dynamically unsta- 
ble mass transfer. Embedded in the CE are the 
dense core of the primordial primary and the pri- 
mordial secondary. However, if the primordial sec- 
ondary is only slightly less massive than the pri- 
mordial primary, the CE may obtain enough en- 
ergy to eject itself even though the orbit of the 
new binary decays slightly. Then, it is difficult for 
the system to fulfill the conditions for SNe la. On 
the contrary, because the orbital period of a sys- 
tem with a lower massive secondary must shrink 
heavily to eject CE, the system may fulfill the con- 
dition for SNe la more likely. Actually, almost all 
of the primordial secondaries leading to SN la have 
a mass less than 4.5 Mq (see Section f5. 2. 2[) . 

As noted by IWang et all (|2009tJ ). WD + He 
star channel may contribute to very young pop- 
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ulation of SNe la, possible all young population 







a CE =1.0 : 


o 




a CE =3.0 ] 








CM 








o 
















CO 






/ ' \ 


o 






/ / V i 










r I ' *V 












o 








i 1 v 










1 1 /\\ : 










i v / \ \ : 


ID 

o 





















8 9 10 

log(«/yr) 

Fig. 7. — The evolution of the birth rates of SNe la 
for a single starburst of 10 11 M Q for different «ce 
(solid lines: «ce = 1-0; dashed lines: «ce = 3.0). 
The thick lines are the results in this paper, while 
the thin lines are the results from the WD + He 



star channel from lWang et al.1 (|2009bh 




Fig. 8. — The evolution of the birth rates of SNe 
la for a single starburst of lO^Mp, for qqe =1.0 
including WD + He star (|Wang et alj2009bl ). WD 
+ MS and WD + R G channel. The points are from 
Totani et al. I (l2008h and the crosses represent error 
bar. 



Using t he results in this paper and in lWang et al 
(2009b), we constructed the evolution of the birth 
rates of SNe la for a single starburst of 10 u Mq for 
«ce = 1.0, shown in Fig. [5] The constructed dis- 
tribution of delay time (DDT) is much similar to 



that d erived from observations bv lMannucci et al 



( 20061 ). except that the peak valu e of young pop 



ulatio n is smaller than that in iMannucci et al 
(|2006h . The WD + He star channel produces 



10 per cent of all SNe la, which is the wea k bi- 
modality as suggested by M annuccil (J2008). In 
the figure, we also showed the observational re- 
sults from iTotani et all (|2008l) . The results in 



th is paper seem to be slightly smaller than those 



in ITotani et al.l (|2008l ). How ever, consider i ng th e 
large error of observation in ITotani et all (|2008t ). 
our results are not inconsistent with observations 
at least. 

5.2. Distribution of Initial Parameters of 
SD Systems for SNe la 

Observationally, some SD systems are pos- 
sible progenitors of SNe la (see the review of 
Parth asarathv et al. I l2007h . Further studies are 
necessary to finally confirm them (from both ob- 
servations and theories). In this section, we will 
present some properties of initial SD systems for 
SNe la , which may help us to search for the po- 
tential progenitors of SNe la in the Galaxy. 

5.2.1. Distribution of Initial Masses of WDs 

In Fig. HI we show the evolution of the distri- 
bution of the initial WD masses in SD systems 
for different «ce- We also show the evolution 
of the upper and lower boundary of the distri- 
bution. From the figure, we can see that most 
of WDs explodes as SNe la in 2 Gyr later af- 
ter starburst, where WD + MS systems are the 
main contributor. acE does not significantly af- 
fect the evolution of the distribution except for 
WD + MS systems with age l ess th an 2 Gyr. As 
shown in Meng. Chen fc Han (l2009h . a high acE, 
i.e. acE = 3.0, leads to the disappearing of the 
WD + MS system from TPAGB channel. So, we 
can not see high-mass WDs for acE = 3.0 when 
age is less than 2 Gyr. For the case of q?ce = 3.0, 
almost all the high-mass WDs are from WD + RG 
channel. There is a clear trend that the range of 
the WD mass decreases with age, while the mean 
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value of WD mass increases. We will discuss this 
in the Section [6l 

5.2.2. Distribution of Initial Secondary Masses 

Fig. [TU] presents the evolution of distributions 
of the initial secondary masses for SNe la with 
different «ce- We also show the evolution of the 
upper and lower boundary of the distribution in 
the figure. It is clear that the distributions of the 
initial orbital period and the evolution of the dis- 
tribution for different acE are similar. The upper 
limit of the secondary mass is 4.5Mq, lower than 
5Mq obtained from binary evolution calculation, 
which indicates that the mass-stripping effect can 
not contribute to the young population of SNe la. 

5.2.3. Distribution of initial orbital periods 

Fig. [TT] presents the evolution of distributions 
of the initial orbital period for SNe la for different 
cxce ■ We also show the evolution of the upper and 
lower boundary of the distributions in the figure. 
It is clear that the distributions of the initial or- 
bital period and the evolution of the distributions 
for different «ce are similar. 
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Fig. 9. — The evolution of the distributions of the 
initial WD masses in SD systems for different oce, 
where a single starburst is assumed. The dotted 
lines shows the evolution of the upper and lower 
boundary of the distributions. 



DISCUSSIONS 



6.1. Comparisons with Previous Studies 



Previous studies show that the minimum mass 
of CO WDs leading to SNe la, 
low 



MSS, may be 



as 



as 



0.70M W for Z = 0.0 2 (|Langer et al 
20001 lHan fc Podsiadlowskilliool . and this min- 
i mum mass even heavily d epends on metallicity 
(jMeng. Chen fc Hani [20091) . In this paper, the 
minimum mass is as low as O.565M0, which is 
much lower than 0.70M Q . This is because that 
the mass-stripping effect and the effect of thermal 
instabl e disk are incorporated in our study. Re- 
cently, IWang. Li fc Hani (j2009l) carried out a se- 
ries of very similar binary evolution calculation to 
this paper and obtained some much similar results, 
especially for the WD + RG channel which are al- 
most the same as those in this paper. In their 
paper, the minimum mass limit is about 0.61 M Q , 
higher than that in this paper. This difference re- 
sults from the mass-stripping effect considered in 
this paper. 

In Fig. [12J we show a comparison of the 
results in this paper with those obtained in 
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Fig. 10. — The evolution of the distributions of 
the initial secondary masses in SD systems for dif- 
ferent acE, where a single starburst is assumed. 
The dotted lines shows the evolution of the upper 
and lower boundary of the distributions. 
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Meng. Chen fc Hani (|2009l ). Our results 
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Fig. 11. — The evolution of the distributions of the 
initial orbital periods in SD systems for different 
ctCE, where a single starburst is assumed. The 
dotted lines shows the evolution of the upper and 
lower boundary of the distributions. 




Fig. 12. — A comparison o f the results in this 
paper with those obtained in iMeng. Chen fc Han 
(2009). Dark contours show the parameter regions 
in the (log P^M^) plane for different WD masses 
that lead to a SN la in thi s paper. The light 



grey c ontours are taken from IMeng. Chen fc Han 
(|2009h . 



much different from those in IMeng. Chen fc Han 
(l2009h. even only f or WD + MS channel. In 
Meng, Chen k, Hani ([2009), their upper limit for 
donor stars of 3.7M Q in the progenitor bina- 
ries is substantially smaller than the limit of 5.0 
Mq obtained here, which arises from the mass- 
strip ping effect. The l ower limit for donor stars 
Meng. Chen fc Hani (|2009h is also much higher 



than those in this paper, which results from the 
thermal instable disk model. The upper limits of 
donor stars for various initial W D mass in this pa- 
per ar e even higher than those in I Wang. Li fc Han 
(j2009h . which arises from the consideration of the 
mass-stripping effect in this paper. 

In addition, WD + RG channel are also con- 
sidered in this paper. In our results, there is no 
WD + RG system wi th period as long as 10 2 - 10 3 
day s as indicated i n iLi fc van den Heuvell (Il997h 
andlNomoto et all (|l999l 120031 ) (jWang. Li fc Han 
(|2009h obtained a similar result to ours). This is 
because if the period of an initial binary system is 
too long, the mass transfer rate between WD and 
donor star is too high and optically thick wind will 
occur and take much hydrogen-rich material away 
from the binary system. The donor star then has 
no enough material to accumulate on to WD. 

6.2. Other Possible Channels and Mecha- 
nisms for SNe la 

Fig. [6] shows Galactic birth rates of SNe la for 
the SD channel. The simulations give a Galac- 
tic birth rate of 2.55 — 2.9 x 10~ 3 yr _1 (including 
WD + He star channel), which is almost consistent 
with but seems to be slightly lower than the birth 
rate inferred observationally. Our results are als o 
slightly smaller than those in lTotani et al.l (|2008h . 
In Fig. [13l we compared our results (including 
WD + He star chann el) with the observation from 
Sullivan et all [2 006). We obtained a similar trend 
to that in Sullivan et al.l ( 2006 ). i.e. the birth 
rate of SNe la is proportional to SFR. However, 
the comparison also shows a slightly lower result 
than observation, while the example normalized 
by Galactic birth rates of SNe l a may be well 



consistent with the observation in Sullivan et al 



(2006J). Therefore, there may be other channels or 
mechanisms contributing to SNe la. 

Although, the WD + RG system is included 
here, it is still should be studied carefully in 
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the future, since this channel may explain some 
SNe la with long dealy time. In addition, 
RS Oph and T CrB (recurrent nova, belong- 
ing to WD + RG) are both su ggested to be 



candid at es of SNe la pr o genitor (Hachisu et al 



1999b] lHachisu fc Katol l2006bt lHachisu et al 



20071 ). while have a muc h longer orbital period 



than those in this paper ( Lines, Lines & McFaul 
1988[lDobrzvcka fe KenvorJl994UDobrzvcka et al.l 
1996). Perhaps, Symbiotic binaries, consisting 



of WD and (super)giant, not filling its Roche 
lobe and havi ng stellar wind sho uld be consid- 
ered carefully ( Branch et al. 19951 ). Additionally, 
the progenitor of some SNe la (e.g. SN 2006X 
and S N 2007on) are possible WD + RG sys - 
tems (|Patat et al.l 120071 IVoss fe Nelemansl 120081) . 
although WD + MS channel may also be the 
candidates for the progenitor of the SNe l a 
(|Hachisu et al.ll2008l : iMeng. Yang fe Geng|[2009ah . 



An alt ernative is the double- d egenerate (DP) 
channel (|lben fe Tutukovl Il984 IWhelan fe Iben 



1987) , although it is theoreti cally less favored 
(Hill ebrandt fc Niemeven HoOO). In this chan- 



nel, two CO WDs with a total mass larger 
than the Chandrasekhar mass limit may coa- 
lesce and explode as a SN la. The birth rate 
from this chann el is comp ar able to the obser- 
vational v a lue (Han 1998 : Yungelson fc Liviol 

19981 120001 : iTutukov fc Yungelsonl I2002I ) . and SN 



2003fg and S N 2005hi likely res u lted from the 
DP channel dHowell et al.l |2006f iBranchl 120061: 
iQuimbv. Hoflich fc Wheeler! 120071) . Observation- 
ally, a large amount of PP s ystems are discov- 
ered (jNapiwotzki et al.l l2004h . but only KPP 
1930+2752 is a possible pr ogenitor candidat e 
for a SN la via PP channel (jGeier et al.ll2007l) . 
The total mass of KPP 1930+2752 (~ 1.52M Q ) 
exceeds the Chandrasekhar mass limit and the 
time scale of coalescence is about 200 Myr esti- 
mated from orbital shr inkage caused by gravita- 
tional wave radiation ( Geier et"aD 2007|) . How- 
ever, lErgma. Fedorova fc Yungelson! ( 2001 ) ar- 
gued that, from detailed binary evolution cal- 
culation, the final mass of KPP 1930+2752 is 
smaller than the Chandrasekhar mass limit due 
to a large amount of mass loss during evolu- 
tion. In addition, earlier numerical simulations 
showed that the most probable fate of the coa- 
lescence is an accretion-induced collapse and, fi- 
nally, neutron star formation (see the review by 



Hillebrandt fc Niemever 2000f) . A definitive con- 
clusion for PP model is thus premature at present, 
and further studies are needed. 

In this paper, we have not considered the influ- 
ence of rotati on on the H-accreting WPs. The cal - 
culations by Yoon. Langer fc Scheithauerl (|2004l ) 
have shown that the He-shell burning is much 
less violent when rotation is considered. Then, 
the He-accretion efficiency may be significantly in- 
crease Meanwhile, the maximum stable mass of 
a rotating WP may exceed the Chandrasekhar 
mass li mit (i.e. the super-Chandrasekharh mass 
model: Ucnishi et al. 2003; Yoon fc Langer 2005; 
Chen fc Li 20091 ). However, we only focus on the 
standard Chandrasekhar mass explosions of the 
accreting WPs in this paper. 

and Wickramasinghe fc Ferrariol 



Liebert et al 



20031 .1200, 

((2005J) found that about 10% of WPs have mag 
netic fields higher than IMG. The mean mass of 
these WPs is 0.93 M©, compared to mean mass 
of all WPs which is 0.56 Mq (see the review by 



Parthasarathv et alJ 120071 for details). Thus, the 
magnetic WPs are more likely to reach the Chan- 
drasekhar mass limit by accretion. Meanwhile, 
the magnetic field may also affect some properties 
of WP+MS systems, e.g. the mass transfer rate, 
the critical accretion rate, the thermonuclear re- 
action rate etc, leading to a different birth rate of 
SNe la. 

6.3. C/O Ratio: the Origin of the Lumi- 
nosity Scatter of SNe la? 

It is widely known that there exists a scatter 
of the maximum luminosity of SNe la and the 
scatter are affected by their environment. The 
most luminous SNe la always occur in spiral galax- 
ies, while both spiral and elliptical galaxies are 
hosts for dimmer SNe la, which lead to a dim- 
mer mean p eak brightness in ell iptical than in spi- 
ral galaxies (jHamuv et al.lll996l ). In addition, the 
mean peak brightness of SNe la in a galaxy has 
less variatio n in the outer regions than in the in - 
ner regions (|Wang et al J ll997b iRiess et al.lfl999h 



To ex plain these phenomena, iNomoto et al.l (| 199 



20031) suggested that the ratio of carbon to oxy- 
gen (C/O) of an initial CO WP is the origin of 
the luminosity scatter of SNe la. The ratio is 
a function of WP mass, i.e. a massive CO WP 
leads to a lower C/O ratio, and thus a lower 
amount of 56 Ni synthesized in the thermonuclear 
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explosion (jNomoto et al. I1999L 12003^ which re- 
sults in a lower luminosity of SNe l a iArnettl 1982 : 
Arnett. Branch fc Wheelen Il985l : iBranchl Il992¥ . 
In Fig. |H1 we showed the evolution of the range of 
WD mass with time and found that the range de- 
creases, while the mean mass of CO WD increases 
with time. If the C/O is the origin of the lumi- 
nosity scatter of SNe la, then our re sults may well 
explai n the phenomenon found by lHamuy et al. 
( 19961 ). In a spiral galaxy, because of existing con- 
tinuous star formation, and then CO WD with 
various age and initial mass, the spiral galaxy may 
be host for luminous and dim SNe la. However, In 
an elliptical galaxy, no star formation and no CO 
WD with young age, only massive CO WD may 
contribute to SNe la, and then a dimmer SNe la 
is expected. The different variation of mean peak 
brightness of SNe la between inner and outer re- 
gion in a galaxy can also be well explained by the 
different WD mass res ulting from different metal - 
licity (see details in Mens. Chen fc Han 2009). 
So, ou r results and those in iMeng. Chen fc Han 
(|2009l ) uphold the C/O as the origin of the lumi- 
nosity scatter of SNe la. Our results in this pa- 
per provide a method to check whether or not the 
C/O is the origin of the luminosity scatter of SNe 
la. If the C/O is the origin, the elliptical galaxies 
with age less than 2 Gyr should be hosts of SNe la 
with various luminosity, including very luminous 
SNe la, while older galaxies are only the host of 
dimmer SNe la. In other words, age might be the 
most important factor to determine the luminosity 
of SNe la, and dimmer SNe la would have a wide 
age distribution. Interestingly, m any observations 
have o btained sim i lar conclusions (iGallagher et al " 
2008; iNeilll 120091: iHowell et all |2009b|). Consid- 



ering the mass of secondary decreases with age 
(see Fig. I10[) . a progenitor with more massive 
secondary shou ld be correlated wit h more lumi- 
nous explosion (IHowell et ah 2009b ). In addition, 
IHowell et al.l ( 2009bl ) also noticed that low-stretch 
(dimmer) SNe la be evenly distributed between 
host age older than 5 x 10 s yr, which is consistent 
with our results qualitatively. Our results may 
also indicate that the mean maximum luminosity 
of SNe la should increase with redshift. Interest- 



i ng, b oth IHowell et al.l (|2007| ) and ISullivan et al.l 
(2009) bund that the mean 'stretch factor' in- 



creases with redshiftQ. 

6.4. Young Population of SNe la 

Observationally, there exists very young SNe la 
(younger than 10 s yr). To explain thes e SNe la in 



Hachisu et all (|2008l ) de- 



the frame of SD model 
signed a WD + MS channel with mass-stripping 
effect. However, the detailed BPS study in this 
paper showed that the WD + MS channel with 
mass-stripping effect totally can not produce SNe 
la with very short delay time. The WD + He 
star channel can w ell produce the young SNe la 
(IWang et al.ll2009al Tbh . The WD + MS, WD + RG 
and WD + He star channel can produce a weak 
bimodality of DDT, which might mean that there 
still exist other channel contributing to the young 
SNe la. Maybe, a symbiotic star with aspheri- 



cal st ellar wind is a possible candidate (jLii et al 
20091) . 



7. SUMMARY AND CONCLUSIONS 



Inc orporating mass-stripping effect in lHachisu et al.l 
(|2008l ) and the effect of the instability of accre- 
tion disk on the evolution of WD binaries into 
Eggleton's stellar evol ution code, and includ ing 
also the prescription of Hachisu et al.l ( 1999a ) for 
mass-accretion of CO WD, We performed binary 
evolution calculations for more than 1600 close 
WD binaries. Adopting the results obtained here, 
we carried out a series BPS sutdy and calculated 
the birth rate of SNe la. We summarize the basic 
results as following. 

1. The detailed binary evolution calcula- 
tion results further confirmed the suggestion in 
Hachisu et al.l (|2008l ) that the mass-stripping ef- 



fect may attenuate between WD and the compan- 
ion, and avoid the formation of a CE, and then 
increase the donor mass leading to SNe la. For 
a reasonable strength of the mass-stripping effect 
in physics, i.e. c\ — 1.5, a companion as massive 
as 5 Mq can produce an SNe la. However, the 
detailed BPS study show that the upper limit of 
the companion could be ~ 4.5M Q . 

2. The detailed binary evolution calculation re- 
sults confirm that the disk instability could sub- 
stantially increase the mas- accumulation efficiency 



4 Stretch factor is an indicator of the luminosity of SNe la. 
A luminous SN la usually has a higher stretch factor. 
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for accreting WD, and extend the mass of donor 
star p rod ucing SNe la to ~ LMq (see also lXu fc Li 
l2009l and lWang. Li fc Hanll2009l ). 

3. CO WDs leading to SNe la may have a mass 
as low as O.565M . 

4. The Galactic birth rate from the SD chan- 
nel in this paper is 2.25 — 2.6 x 10~ 3 yr" 1 , based 
on acE- If the WD + He star channel is also in- 
cluded, the Galactic birth rate from the SD chan- 
nel increases to 2.55 — 2.9 x 1CP 3 yr _1 , only slightly 
smaller than that derived from observation. Then, 
there should be other channels or mechanisms con- 
tributing to SNe la. 

5. The DDT from the WD + MS, WD + RG 
and WD + He star channels may be a weak bi- 
modality. The WD + MS produces most of SNe 
la with delay time between 0.15 — 2.5 Gyr. The 
WD + RG channel may mainly contribute to SNe 
la with delay time longer than 6 Gyc. 

6. The mass-stripping effect may not contribute 
to very young population of SNe la. The WD + 
He star channel could contribute most of young 
SNe la, even though not all. 

7. Our BPS results may uphold the C/O as 
the origin of the luminosity scatter of SNe la and 
provide a method to verify it. 

8. We show the evolution of the initial param- 
eters of SD systems for SNe la with time, which 
is helpful for searching for possible candidates of 
SNe la progenitor systems. 
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Fig. 13. — SNe la rate per unit stellar mass 
as a function of the SFR per unit stellar mass 
fo r various ap y . The WD + He star channel 
in IWang et all (|2009bh is included. The dotted 



line represents an example normalized by Galac- 
tic b irth rates of SNe la. The data are obtained 
from Sullivan et al. (2006) and the cross represent 



the error bar. 
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